Poisson noise induced switching in driven micromechanical resonators 
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We study Poisson-noise induced switching between coexisting vibrational states in driven non- 
linear micromechanical resonators. In contrast to Gaussian noise induced switching, the measured 
logarithm of the switching rate is proportional not to the reciprocal noise intensity, but to its log- 
arithm, for fixed pulse area. We also find that the switching rate logarithm varies as a square root 
of the distance to the bifurcation point, instead of the conventional scaling with exponent 3/2. 
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I. INTRODUCTION 

Non-Gaussian fluctuations in mesoscopic devices and 
in photon statistics have proved to be a powerful probe of 
the underlying physics of such systems, yielding informa- 
tion that cannot be deduced from the averaged value and 
the varianceii Experimentally, detecting a deviation from 
the Gaussian distribution is a challenging task. Direct 
measurements of the third and higher moments of fluctu- 
ations in mesoscopic systems require highly sophisticated 
techniques 1^"— In another approach, following theoretical 
proposals r'"— non- Gaussian statistics was detected using 
noise-activated switching out of a metastable state in a 
Josephson junctionJ^iii This approach provides direct 
access to the tail of the noise distribution, as a relatively 
large noise outburst is required to drive the Josephson 
junction out of the effective potential well. In the mea- 
surements, deviations from the Gaussian statistics were 
found from the asymmetry of the switching rate with re- 
spect to the current polarity. 

Noise-induced switching has attracted much attention 
recently in the context of nonequilibrium systems, where 
the notion of escape from a potential well does not apply. 
In particular, a number of experiments were done on peri- 
odically driven systems, with the metastable states corre- 
sponding to forced vibrations. The systems ranged from 
trapped electronsi^ to Josephson junctionsr^ mechani- 
cal resonators r^ii^ and trapped atomsJ^ It was found 
that, for Gaussian noise, the switching rate W displays 
universal features. Its dependence on the noise intensity 
Dg is of the activation form, W oc exp(— Q), with the 
switching exponent Q = Rg/Dq, where Rq is the ef- 
fective activation barrier. Furthermore, Rg displays the 
predictedii scaling behavior close to bifurcation points 
where the number of stable states changes fi^"— Rg oc 77^, 
where rj is the distance to the bifurcation point. 

In this paper we use a driven micromechanical oscilla- 
tor to study the features of Poisson-noise induced switch- 
ing. In contrast to Gaussian noise, Poisson noise is char- 
acterized not just by its intensity Dp and correlation 
time, but also by the pulse rate lyp. We limit ourselves 
to the most interesting case of short pulses and explore 



the dependence of the switching rate on the pulse area 
and vp as well as the distance to the bifurcation point 
rj. Even though the stable states of this nonequilibrium 
system are not separated by a potential barrier, our find- 
ings are consistent with the recently predicted universal 
behavior of Poisson-noise induced switching near a bi- 
furcation point, which differs qualitatively from that for 
Gaussian noissi^. 
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FIG. 1. (a) Schematic of the experimental setup. The noise 
force is generated by applying voltage Vdc + 'Kioiso to the top 
sidegate, with |Vdc| 3> |Vnoisc|. (b) Vibration amplitude vs. 
driving frequency; the vibration ampUtude is scaled by its 
value at the critical point where hysteresis emerges. The 
curves correspond to the driving amplitudes of 0.2^c, Ac, and 
4:Ac, where Ac is the critical driving amplitude for the onset 
of hysteresis. The up-arrow indicates the position of the bi- 
furcation point cji. We fit the linear and the critical responses 
to extract the parameters of the oscillator, (c) Schematic of 
the generation of Poisson pulses. Each output pulse contains 
> 400 rf cycles. For clarity, only 2 rf cycles are shown. 



II. EXPERIMENTAL SETUP 

Our oscillator is a double-clamped polysilicon beam 
that is 100 fim long and 1.2 ^m by 1.5 ^m in cross section. 
By passing an ac current through the beam in a perpen- 
dicular magnetic field of 5 T, we excited motion of the 
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beam in its in-plane fundamental mode [Fig. 1(a)]. This 
motion generates an electromotive force that changes the 
transmitted ac power. All measurements were performed 
at 4 K and < 10~^ torr. The mode can be modeled by a 
Duffing oscillatorilii^i^, 

+ 2r J + ujlq + ^q"" =AcoSLOpt + fit), (1) 

where q is the normalized beam displacement, F = 96 
rad s~^ is the damping coefficient, wq = 7,133,339 
rad s~^ is the resonant frequency, 7 = 2.3 x 10^ s~^ 
is the coefficient of the cubic nonlinearity, A and ajp are, 
respectively, the amplitude and frequency of the external 
driving force, and f{t) is the noise. 

In the absence of noise, for small A the device responds 
as a weakly damped harmonic oscillator [Fig. [1] (b)]. As 
A increases, the response curve bends toward higher fre- 
quencies because of spring hardening. For A exceeding 
the critical value A^ the oscillator becomes bistable, with 
two coexisting states of forced vibrations in the range 
LUi < LUp < The frequencies W1.2 depend on A. They 
correspond to saddle-node bifurcations and arc indicated 
by arrows in Fig. [1] (b). 

Noise can lead to switching between the vibrational 
states. In the parameter range of our operation, thermal 
noise was too weak for switching to be observed. We 
generated noise externally, by applying random electro- 
static forces between one of the sidegates and the beam. 
To create Gaussian noise, f{t) = fait) in Eq. ([T]), the 
Johnson noise from a 50 51 resistor at room temperature 
was amplified and directly applied to the sidcgate. 

Generating Poisson noise involves extra circuitry 
[Fig. [1] (c)]. Since the oscillator is most perceptive to 
noise at frequencies close to cjq, we used noise that con- 
sisted of pulses at frequency ujf with duration tp ^ ojp^. 
In the rotating frame, as we will explain later, such pulses 
looked rectangular. To create them, the Gaussian noise 
voltage was connected to the trigger input of a pulse 
generator. Whenever the input voltage, in a rare large 
outburst, exceeded a threshold value, the generator pro- 
duced a square pulse. The pulse width tp = 400 ^s 
was much less than the reciprocal pulse rate Vp^, which 
ranged from 30 ms to 200 ms. We verified that the inter- 
pulse intervals obeyed the Poisson statistics. The pulse 
signal controlled an rf switch that only turned on when a 
pulse was present, generating an electrostatic force with 
the same phase as the regular magnetomotive force. The 
output of the switch is described by the random force in 
Eq. (dl of the form 

f{t) = fp{t) = fp{t)cosu:pt, (2) 

where Jp{t) is a random pulse train, with the area of each 
pulse gp controlled by a variable attenuator; for chosen 
tp there were « 454 rf cycles in each pulse. The pulses 
represent Poisson noise with intensity Dp — vpg'^/2. We 
checked that a single pulse never led to switching. 



III. THEORETICAL BACKGROUND 

To study the oscillator dynamics, we go to the rotat- 
ing frame and introduce slow variables X and Y , q{t) = 

Cresi^ cos ujpt + Y sin ujpt), q{t) = —ujpCres(X sinujpt — 

1/2 

Y cosujpt), where Crcs = {SiOpSuj/Sj) and 6uj ~ 
LUp ^ loq is the frequency detuning. Functions X and 

Y are the scaled in-phase and quadrature components of 
q{t). In the rotating wave approximation (RWA) Poisson 
noise ([2]) is demodulated and becomes a random force fp 
that drives the quadrature component Y. In dimension- 
less "slow" time r = tSuj = t{ujp — wq) Eq. ([T]) becomes 

A = Ax, Y = KY+~fp{T), 

Kx = -n-^x + OyG, Ky = -n-^Y - dxQ, 

g{X,Y) = ^{X^ + Y'^ -If - p^/^X. (3) 

Here, dot indicates differentiation with respect to r, and 
the noise fpir) = f p{t) / 2uj pC-^as^uj . The two dimension- 
less parameters (3 and are the scaled intensity of the 
modulation and the scaled detuning of the modulation 
frequency from the oscillator eigenfrequency, 

/3 = 37^V32w|((5w)^ VL-'^^T/5oj. (4) 

Functions Kx , Ky do not explicitly depend on time in 
the rotating frame. The stationary solutions of equations 
(|3|) in the absence of noise, which are given by equations 
Kx = Ky = 0, determine the stationary vibrational 
states of the oscillator. In the region of bistability there 
are three stationary states, two of them are stable and 
the third is the saddle-type solution of Eq. ([3]) in the 
absence of noise. 

For chosen pulse duration tp <^ F^^, the noise fpir) 
is effectively (5-correlated. From Eqs. ([T]) and each 
pulse of /p(t) shifts Y{t) by gp — {2ujpCycs)^^gp- The 
dimcnsionlcss pulse rate is Dp — vp/Suj and the dinicn- 
sionless intensity is Dp = vpgp/SujpC'^^gSoj. 

If the system is initially at a stable state, for small 
vp/T and gp it usually has time to relax back to this 
state between the pulses. However, occasionally there 
may happen a sequence of frequent pulses that will push 
the system far from the occupied state so as to cause 
switching to a different state. This switching mecha- 
nism is qualitatively different from the case of Gaussian 
noise where individual noise pulses are too frequent to 
be resolved by the system and the noise outburst has to 
overcome the deterministic force that drives the system 
toward the stable state. 

It is most interesting to study Poisson-noise induced 
switching close to the saddle-node bifurcation point 
where, on the one hand, the switching rate W becomes 
larger and, on the other hand, the system dynamics is 
expected to display universal features related to the crit- 
ical slowing down. The quadrature Y(t) corresponds to 
the soft mode near a bifurcation point of the Duffing 
oscillator. The in-phase component X{t) follows Y{t) 
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adiabatically, so thaljii 

X{T)-XB^aB[Y{T)-YB].. 



as = f^(2r|-l), (5) 



where Xb and Yb arc the values of X and Y at the 
bifurcation point and = X^ + . 

Extending the analysisii to the Poisson noise, one 
obtains equation of motion for the soft mode as Y = 
-dU/dY + fp{T). The potential U has the form U{Y) = 
ciy^ + C2i]y, where y = Y — Yb is the distance to the 
bifurcational value of Y. The parameters ci 2 in the ex- 
pression for U arc knownjiS 77 is the distance to the bifur- 
cation point in the space of the driving force parameters 
{A, lof)', If ciC2?7 < 0, potential U (Y) has a minimum and 
a maximum at l^in and ymax, respectively. The mini- 
mum corresponds to the stable state of forced vibrations. 
The extrema of U{y) merge at the bifurcation point 77 = 
and disappear for ciC2?7 > 0. This maps the problem of 
the oscillator dynamics near a bifurcation point onto the 
problem of an overdamped particle with one dynamical 
variable in a static potential. 

Using the results for the switching rate W oc exp(— Q) 
near a bifurcation point in a static potential;^, we obtain 
for the switching exponent Q = Qp 

Qp « (2,71/7,9) log(K77/gi.), v = I3b- P- (6) 



Here 



9 ■ 



2ujfCb 



9P, 



5u}R 



-Vp, 



(7) 



with h = /3^'^((5wB)2(3r|-2)/2r2. Parameters /3b, Sujb, 
and Cb are the bifurcational values of the scaled modula- 
tion intensity /3 given by Eq. (|4]), the frequency detuning 
5u), and the scaling factor Cics, respectively; parameter k 
in Eq. ([5]) is 1 . Equation ([5]) applies provided the num- 
ber of noise pulses rigw needed to shift the system from 
l^in to l^iax is n-sw 3* 1, {^p / 5uj)t-c^ whcrc the dimension- 
less relaxation time is — 1/U" {Ynun)- We checked that 
this condition was met in the experiment. Intuitively, 
Qp can be found by maximizing the probability for the 
Poisson distribution to have Usv, ^ (i^iax — i^min)/5P <x. 
T^l^ I g pulses in time Tj. and taking the logarithm of this 
probability!^ 

Equation ([6]) should be contrasted with the expression 
for the switching exponent Q = Qq in the case of Gaus- 
sian noise, j(t) = /g(0 in Eq. ([1]). For white noise, 
(/G(0/G(t')) = 2DG5{t - t'), near a bifurcation pointii 



3/2 



' (jjj^CgSuJB 



(8) 



This expression displays a different dependence on rj and 
the noise parameters than Eq. ([5]). 



IV. RESULTS AND DISCUSSION 

In the experiment, the driving amplitude A was ~ 4:Ac, 
so that the hysteresis loop in Fig. [T] (b) was compara- 
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FIG. 2. (a) Transitions from the low-amplitude state to the 
high-amplitude state at Acj,, = 9.4 rad s~^ in the presence 
of Poisson noise. The amplitude is normalized to the vibra- 
tion amplitude at the critical point. The small spikes are the 
response of the oscillator to pulses that do not lead to tran- 
sitions. The response during reset to the low-amplitude state 
is not shown, (b) The histogram of the residence times. The 
linear fit gives the mean residence time. 



tively large while the vibrations remained almost sinu- 
soidal. We measured the switching rate near the bifur- 
cation point uiB at driving frequency 



OJp — LUlB 
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the distance to the bifurcation point was along the fre- 
quency axis, with Aw,, cx 77. 

The oscillator was initially prepared in the low- 
amplitude state. Then the noise was turned on, and we 
monitored the oscillator amplitude. Figure [5] (a) exhibits 
a time series of the measured amplitude in the presence 
of Poisson noise. After switching to the high amplitude 
state, we reset the oscillator to the low-amplitude state 
by a procedure described in Ref. [l^. Figure [5] (b) shows 
the histogram of the residence time, demonstrating an 
exponential decay of the state population. The linear fit 
determines the switching rate W . A similar behavior is 
observed for switching due to Gaussian noise>i^. The mea- 
sured switching rates for Poisson noise were close to the 
theoretical values that took the prefactor into account^. 

The dependence of the switching rate W on the noise 
parameters and the distance to the bifurcation point is 
presented in Figs. 3-5. Figure[3](a) shows that, for Gaus- 
sian noise, — logM^ linearly depends on the reciprocal 
noise intensity I/Dq, in agreement with Eq. ([5]). The 
Poisson noise intensity Dp = vpgp'^ /2 depends on both 
the pulse rate vp and area gp. In Fig.[3](b) we show that 
the increase of — log with \/vp for fixed is strongly 
sub-hnear. Replacing l/up with log(l/^'p) gives a good 
linear fit, as seen from the inset in Fig.[3](b). This is com- 
pletely different from the 1 /Dq dependence for Gaussian 
noise and is consistent with Eq. ([6]). 

Next, the dependence of the switching rate on the pulse 
area gp was studied by changing the pulse height, while 
keeping the pulse width tp constant. In Fig. 2] (a) we 
show — logW^ as function of log(l/i^p) for different gp. 
Fitting the data by straight lines, we found the slope 
Sp ~ \d\ogW/d\ogi'p\ « \dQ-p/d\ogv\ as function of 
gp. Figured] (b) shows a log-log plot of this slope. The 
measured exponent of Sp oc gp'^ is a = 1.12 ± 0.14, 
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FIG. 3. (a) The logarithm of the switching rate scaled by 
the oscillator decay rate — log(W/r) for Gaussian noise as a 
function of the noise intensity Da for Atj,, — 3.14 rad s~^; 
Dg is in the units of s~^. The solid line is a linear fit. (b) 
The dependence of — log(VK/r) on the reciprocal scaled pulse 
rate V/up at Ao;,, = 9.4 rad s~^ for Poisson-noise induced 
switching; the pulse area gp is constant. Inset: the same 
data plotted vs. \og{T /vp). 

consistent with the theoretical value of 1 in Eq. ^ [the 
factor log gp in Qp, Eq. ([5]), could not be measured in 
the studied range of gp]. 

Because the pulses of fp{t) "push" the quadrature Y{t) 
in one direction, near a bifurcation point, where the mo- 
tion in the rotating frame is ovcrdampcd, switching is 
possible only where this is the direction from Ymin to 
^max- For gp > 0, Eq. ^ describes switching near wis 
in Fig. [1] Switching near the second (larger-frequency, in 
our case) bifurcation point UI2B occurs if < 0. This 
was tested in the experiment by changing sign of fp{t) 
in Eq. The switching exponent near Ltj2B is given by 
Eq. ^ with gp — > \gp\ and rj —5- 
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FIG. 4. (a)The logarithm of the scaled switching rate 
— log{W/T) vs. the logarithm of the reciprocal pulse rate 
log(r/:^p) at Aoj^ = 9.4 rad s^^. The solid lines are a linear 
fit for a fixed gp. From top to bottom, the values of gp were 
3, 3.5, 4, 4.5, 5, 6, in scaled units. The error bars are smaller 
than the markers, (b) Logarithm of the slopes from (a) as 
a function of \og{l /gp). The line is a linear fit, yielding the 
exponent a = 1.12 ± 0.14. 

Last but not least, we studied the dependence of the 
switching exponent on the distance to the bifurcation 
point. For Gaussian noise, we generated a number of 
plots similar to Fig. [3] (a), each at a different Auj,f. We 
could then extract the activation barrier Rq = QgDg 
[the switching rate W cx exp(— i?^/!)^)] from the slope 
of Qg vs Dq^. Figure [5] (a) shows logi?G vs logAcj^. 
A power law fit gives a critical exponent of 1.31 ± 0.20 
in agreement with the theoretical predictions of 3/2 in 
Eq. ([5]) and with other expcriment a^'^i^^ . 



For Poisson pulses, the slopes Sp were extracted from 
plots of —\ogW vs. \og{l/vp) [similar to the inset of 
Fig. El (b)] at different Aw,, and plotted as Fig. [S] (b). 
Fitting the results in Fig. [5] (b) with a power law 
Sp K, \dQp / d\ogvp\ cx rf' gives a critical exponent of 
^ — 0.61 ± 0.08, consistent with the theoretical predic- 
tion of 1/2 [Eq. ©]. 

Apart from the different exponent in the dependence 
on Ao;,, compared to the Gaussian noise, the dependence 
of Qp on /S.ujri for Poisson pulses contains an extra log 
factor in Eq. (jB]). The dynamical range of the parame- 
ters in our experiment did not allow us to single out this 

1 /2 

factor on top of the Aw,, directly by measuring log W 
for different Aw,, while staying close to the bifurcation 
point. However, we made use of the fact that vp is only 
present inside the log factor in Qp. Therefore the slope 
Sp should be independent of log Aw,,. By varying the 

1 /2 

Poisson pulses via vp, the factor Aw,,' could be isolated 
using the slope of the activation exponent in the inset of 
Fig. m (b) as function of log vp for different modulation 
frequency. 

The parameters of the oscillator used in the experi- 
ment were chosen to be close to the bifurcation values. 
However, they were not close enough to make the oscilla- 
tor motion in the rotating frame fully controlled by one 
dynamical variable Y{t)^ i.e., the assumption that the 
variable X{t) follows Y{t) adiabatically was not strictly 
valid. In Appendix we show that the full theory of the 
escape rate in the studied range gives results very close to 
the asymptotic expression This is because X{t) be- 
haves nonadiabatically only very close to the metastable 
state of the oscillator, but the contribution of this area 
to the switching exponent Qp is small. 
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FIG. 5. The activation barrier Rq for Gaussian noise induced 
switching (a) and the slope Sp = j d log W/d log i^pj for Pois- 
son noise induced switching (b) as functions of the distance 
to the bifurcation point Acj,,. The fitted critical exponents 
are 1.31 ±0.20 for Gaussian noise in panel (a) and 0.61 ±0.08 
for Poisson noise in panel (b). 



V. CONCLUSIONS 

Our results demonstrate that driven nonlinear sys- 
tems can serve as bifurcational detectors of non-Gaussian 
noise, based on the sensitivity to the noise statistics of 
the dependence of log W on the distance to the bifurca- 
tion point. The detectors can be implemented with nano- 
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mechanical systems, rf Josephson junctions, or other sys- 
tems. A stringent condition is that within the relaxation 
time of the system the mean number of pulses needs to 
remain small. As a rough estimate, to measure fluctu- 
ations in mesoscopic systems with current ^ 1 pA, the 
relaxation time of the detector needs to be ^ 100 ns. A 
nanomechanical detector will therefore require a state-of- 
the-art GHz resonator with a quality factor > 100. 

At the same time, non- Gaussian noise plays an in- 
creasingly important role in mechanical devices as they 
are getting smaller. Besides the discreteness of the 
fluctuating electric charge coupled to the mechanical 
motion^, it can come from attachment or detachment 
of molecules22r^ or random spin flip a^^'^^ . The effects of 
non-Gaussian noise on nano-mechanical systems warrant 
further investigation. 

In summary, we have measured the rate W of Poisson 
noise induced switching in a nonlinear resonator near a 
bifurcation point where the number of stable vibrational 
states changes. The observed dependence of log W on 
pulse rate vp and area gp, as well as on the distance to 
the bifurcation point 77 is consistent with the predicted 
scalingi^, which is strongly different from the scaling for 
Gaussian-noise induced switching. 

This work was supported by the NSF through grants 
No. CMMI-0856374 and 0900666. H.B.C. is supported 
by Shun King Solid State Clusters Lab and HKUST 
600312 from the Research Grants Council of Hong Kong 
SAR. 



Appendix A: THE MOST PROBABLE 
SWITCHING PATH 

The random trajectories followed by our periodically 
driven micromechanical resonator in fluctuation-induced 
switching are concentrated around a specific path in its 
phase space {X,Y), which for Gaussian noise has been 
now seen in experiment i21 This path is commonly known 
as the most probable switching path (MPSP). Even 
though the mechanism for switching induced by Pois- 
son noise is qualitatively different from Gaussian noise 
as discussed in the main text, the notion of MPSP still 
applies. 

In general, the MPSP follows complex patterns in mul- 
tivariable systems. For example, in a micromechanical 
resonator driven into parametric resonance, if the mo- 
tion in the rotating frame is underdamped the MPSP 
has been shown to consist of spirals around the two at- 
tractors on the (A, Y) plane. As the system approaches 
the bifurcation point, the motion near the corresponding 
stable state in the rotating frame undergoes a crossover 
from 2D to effectively ID and becomes overdamped. In 
a switching event, the system most likely moves along a 
straight line on the (A, Y) plane from the initially occu- 
pied state to the saddle point ; this motion is controlled 
by a ID overdamped soft mode fluctuating in a ID poten- 
tial well. The scaling of the switching exponent predicted 



in Rcf. [1^ is obtained using the asymptotic theory that 
considers the escape of the system from such potential 
well. 

Here, with the parameters used in the experiment, we 
analyze theoretically the motion of our system in phase 
space. We show that under the experimental conditions, 
the portion of the MPSP that goes from the stable state 
to the saddle point (the "uphill" portion) is largely lin- 
ear. Spiraling only occurs close to the initially occupied 
attractor. For the spiraling to entirely disappear so that 
the motion is truly ID, it is necessary to choose driving 
frequencies even closer to the bifurcation value. Never- 
theless, the spiraling portion of the path contributes only 
a negligible amount to the action in the calculation of the 
switching exponent. As a result, the scaling behavior of 
the switching exponent is found to extend beyond the ID 
regime, in agreement with our experimental findings. 
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FIG. 6. (Color online) (a) The uphill portion of the MPSP 
that is numerically calculated with the parameters Aw,, = 
9.4 rad s~^, i/p = 11.5 Hz, and the dimensionless pulse area 
gp = 3.5. The MPSP is shown by the thick solid line (red) 
that goes from the attractor (dot) to the saddle point (star). 
On this scale, the spiraling of the MPSP is not resolved. The 
thin (black) lines are the portions of the separatrix which 
bounds the basin of attraction of the considered attractor. 
The other attractor (not shown) is far away from the sad- 
dle point, (b) Close-up of the MPSP near the attractor. A 
spiraling behavior can be observed. 

In the rotating frame, the generalized Fokker-Planck 
equation for the probability density p of the system in 
the presence of Poisson noise reads^^ 

p = -9q[Kp(q)] + i>p[p(q - gp) - p(q)], (Al) 

where q = (X, Y) is the position vector in the rotating 
frame, K = {Kx , Ky), and gp = (0,gp). The com- 
ponents Kx.Y are given by Eq. gp and iyp are the 
dimensionless pulse area and pulse rate, respectively. 

We are interested in the quasistationary solution of 
Eq. (jAl[) . which is formed over the relaxation time and 
persists for times small compared to the reciprocal rate 
of switching from the considered metastable state. Re- 
spectively, we seek the solution as p(q) = exp[— s(q)]. 

We assume that on the tail of the distribution s(q) is 
large, but that it is smooth on the scale ^ |gp|- Then, 
if we keep terms (x dqS and assume them large, but dis- 
regard terms with higher derivatives of s, wc can reduce 
Eq. (|Aip to equation 

i/(q,p)=0, i7 = i>p[exp(pgp)-l]+p.K, (A2) 
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FIG. 7. (Color online) The measured values of — log(W^/r) 
(small solid circles) agrees with the asymptotic theory (the 
black line) for typical experimental parameters: Ao;,, = 
9.4 rad s~^ and gp — 3.5 in scaled unit. The numerical calcu- 
lated Qp (red circles) are plotted with the right vertical axis 
that is shifted by 1 compared to the left one. 



where p = d^s. This equation maps the problem of 
noise-induced switching on the problem of Hamiltonian 
dynamics of an auxiliary system. This system has the 
coordinate vector q and the momentum vector p, its 
Hamiltonian function is lf(q, p), whereas function s is 
its action,^ii^ The uphill portion of the MPSP is a hetc- 
roclinic orbit that connects the fixed points (Qq, pa = 0) 
and (q5,P5 = 0), where q^i and q^ correspond to the 
positions of the attractor and the saddle point of the res- 
onator in the rotating frame. 

The switching exponent Qp is given by the integral of 
p along the MPSP, 



/■qs 

Qp = p dq. 



(A3) 



We have numerically calculated the MPSP taking into 



account that close, but not too close to the fixed point 
(q^, 0), action s(q) is quadratic in q — Similar to 

how it is done for Gaussian noise (cf. Ref. [23 and papers 
cited therein), one can find MPSP using the shooting 
method by finding the "right" direction for the Hamilto- 
nian trajectories of the auxiliary system that go from the 
vicinity of (q^i, 0), so that they approach (q^, 0). 

Figure [6] shows the uphill portion of a MPSP for typical 
values of the experimental parameters. The trajectory 
does not show any spiraling at this length scale. A close- 
up of the region close to the attractor [Fig. [H] (b)] reveals 
that spiraling indeed occurs. Therefore the uphill path 
is, straightly speaking, not ID. However, since the spi- 
raling portion is small and lies in the region of small |p|, 
whereas much of the uphill path remains largely straight, 
the asymptotic theory provides a good estimation of the 
switching rate. 

In Fig. [71 we compare the measured — \og{W/r) with 
the theoretical values calculated from the the asymptotic 
theory for ID dynamics [Eq. ([5])] and the numerical values 
that go beyond the ID theory. The black line represents 
the resuhs for -log(VF/r) = Qp - log(C/r), where Qp 
and the prcfactor C are both calculated from the asymp- 
totic ID theory>i^ The good agreement indicates that the 
asymptotic theory provides a good estimate of the switch- 
ing exponent even beyond the strictly ID regime. 

The numerical results for the switching exponent Qp 
[Eq. (jA3[) ] refer to the same parameters. Since the 
numerical method does not yield the prefactor C, it is 
not possible to perform a direct comparison between the 
numerically calculated Qp and the experimentally mea- 
sured — log W/r. In Fig. [71 , the right vertical axis (cor- 
responding to the numerically found Q_p) is shifted from 
the left vertical axis (corresponding to — logVF/F). The 
plot indicates that the slope of Q_p vs. logF/z^ is in good 
agreement with the results from both the asymptotic the- 
ory and the measurement. 
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